Multi-Physics Modeling of Melting-Solidification Characteristics in Laser Powder Bed Fusion Process of 316L Stainless Steel

In the laser powder bed fusion process, the melting-solidification characteristics of 316L stainless steel have a great effect on the workpiece quality. In this paper, a multi-physics model was constructed using the finite volume method (FVM) to simulate the melting-solidification process of a 316L powder bed via laser powder bed fusion. In this physical model, the phase change process, the influence of temperature gradient on surface tension of molten pool, and the influence of recoil pressure caused by the metal vapor on molten pool surface were considered. Using this model, the effects of laser scanning speed, hatch space, and laser power on temperature distribution, keyhole depth, and workpiece quality were studied. This study can be used to guide the optimization of process parameters, which is beneficial to the improvement of workpiece quality.


Introduction
Under the assistance of additive manufacturing (AM) technology, periodic or statistically homogeneous materials can be used to explore a wider range of interesting structures [1].As a part of the AM technology, laser powder bed fusion (LPBF) is widely used in the manufacturing of metallic material.By using high temperature restricted to local areas, this technology can directly manufacture the metallic powder into any part geometries or complex specimens with some inner degrees of freedom with a layer-by-layer fashion [2][3][4][5].However, this results in a large temperature gradient, which has great effects on thermal stress.During LPBF, the process parameters have great effects on heat transfer characteristics and molten pool morphology in the melting-solidification process of the powder bed [6,7].The inappropriate selection of process parameters may also result in undesirable effects, such as thermal cracking, balling, and porosity [8].
The 316L stainless steel (SS316L) has been widely used in navigation, nuclear reactor, aerospace, and medical environments for many years due to its perfect mechanical properties, such as high tensile strength, low coefficient of thermal expansion, perfect resistance to oxidation at high temperature, creep resistance, and fatigue resistance [9][10][11].In recent years, the SS316L powder has been widely used in AM to create complex shaped structures due to the properties of excellent corrosion resistance, low carbon content, hightemperature resistance, weldability, and good flowability of the powder bed, which makes it well-suited for AM.Moreover, research shows that the nearly full densification part of SS316L can be achieved using the LPBF process [12].While SS316L offers numerous advantages in AM, there are also some limitations and challenges associated with its use in AM, such as high thermal conductivity, which makes it challenging to achieve optimal melt pool; residual thermal stresses and distortion, which may potentially cause warping or distortion; and inhomogeneous microstructure and porosity, which may affect the mechanical properties.To obtain high-performance final parts of SS316L manufactured by the LPBF process, current research interest in LPBF primarily focuses on optimizing process parameters, analyzing microstructures, characterizing mechanical properties, exploring post-processing techniques, and employing computational modeling and simulation.This collective research effort aims to advance our understanding and enhance the quality of SS316L parts produced via LPBF [13][14][15][16][17][18].
During the LPBF process, the SS316L powder bed is rapidly melted, solidified, and even vaporized.The flow behavior in the molten pool can be affected by these phenomena.Thus, the thermal distribution and molten pool morphology will be changed.To analyze the relationship between the LPBF process and the mechanical properties of final parts, some experimental and theoretical research has been conducted by other scholars.Greco et al. conducted LPBF experiments of SS316L powder with constant laser energy density.Via these experiments, the effects of laser power, layer thickness, and hatch space on workpieces' roughness, density, and micro-hardness were studied [19].Röttger et al. manufactured some specimens using SS316L with different types of LPBF devices.Meanwhile, their microstructures were observed by means of scanning electron microscopy.It was found that these specimens possessed similar microstructures, but there was different porosity between the microstructures, which makes the specimens have significant differences in mechanical strength [20].Suryawanshi et al. compared the part of SS316L produced using LPBF with the conventionally manufactured austenitic stainless steel.By observing their micro and meso-structures, the property variation was analyzed [21].Şener et al. conducted a series of experiments to study the properties of the samples by changing the process parameters.Via optical profilometry and tensile test, the effects of process parameters on mechanical performance and surface quality were obtained and optimized [22].The above experiments can efficiently analyze the mechanical properties of the final parts processed by LPBF, but the effect of the melting-solidification mechanism of the metallic powder bed is hard to analyze for the characteristics of molten pool, such as large temperature gradient and flow behavior at a small scale.
Considering the molten pool characteristics, the mechanism model is suitable to analyze the melting-solidification process caused by the rapid temperature change in the LPBF process.The developed numerical model of LPBF by other scholars can be used to guide the improvement of the mechanism model for SS316L.Waqar et al. developed an FE model of LPBF to analyze the influence of process parameters on thermal and stress profile [23].Trejos et al. also developed an FE model of LPBF, which takes into account most of the manufacturing parameters and thermophysical properties.Using this model, the thermal history of the LPBF process was calculated.By considering nodes with temperatures higher than the melting point as molten status, the evolution of the molten pool dimensions was analyzed [24].Luo et al. developed a new efficient and accurate FE model for predicting transient material thermal behavior during the LPBF process in partlevel.To reduce the computational cost, the adaptive mesh scheme was used, which makes the simulation speed 12× to 18× faster compared to the traditional simulation scheme [25].As the liquid metal flow behavior was not considered in the above models, it is hard to describe the heat transfer process caused by liquid metal flow.Additionally, due to the wide application of LPBF in industrial production, CAE software (ANSYS version 2019R1) providers also developed solutions to simulate the LPBF process.These solutions can accurately analyze the residual deformations and internal stresses under large temperature gradients [26][27][28], but the inhomogeneity of the powder bed distribution and the liquid metal flow make it hard analyze their effects on the melting-solidification process.To analyze the melting-solidification process of SS316L thoroughly, a comprehensive model coupling heat transfer physics and laminar flow physics is eagerly needed.
In this study, a three-dimensional numerical model of LPBF for SS316L was developed to describe the temperature and velocity field in the melting-solidification process.Via the single and multi-track simulation, the effects of process parameters on temperature distribution, keyhole depth, and workpiece quality were studied.By shedding light on the characteristics of the melting-solidification process, the correlation with the process parameters can be clarified.Based on this study, some phenomena during manufacturing can be explained.Meanwhile, it may contribute to the design of experiments when analyzing the correlation between the process parameters and mechanical properties of 316L stainless steel produced by LPBF.

Computational Models
During the LPBF process, complex physical phenomena such as heat transfer, phase change, and fluid flow of liquid metal materials are involved, as shown in Figure 1 [29].By describing the physical phenomena through the corresponding numerical model, the LPBF process can be quantitatively analyzed, which is of great significance to guide the optimization of workpiece quality.
metal flow make it hard analyze their effects on the melting-solidification process.To analyze the melting-solidification process of SS316L thoroughly, a comprehensive model coupling heat transfer physics and laminar flow physics is eagerly needed.
In this study, a three-dimensional numerical model of LPBF for SS316L was developed to describe the temperature and velocity field in the melting-solidification process.Via the single and multi-track simulation, the effects of process parameters on temperature distribution, keyhole depth, and workpiece quality were studied.By shedding light on the characteristics of the melting-solidification process, the correlation with the process parameters can be clarified.Based on this study, some phenomena during manufacturing can be explained.Meanwhile, it may contribute to the design of experiments when analyzing the correlation between the process parameters and mechanical properties of 316L stainless steel produced by LPBF.

Computational Models
During the LPBF process, complex physical phenomena such as heat transfer, phase change, and fluid flow of liquid metal materials are involved, as shown in Figure 1 [29].By describing the physical phenomena through the corresponding numerical model, the LPBF process can be quantitatively analyzed, which is of great significance to guide the optimization of workpiece quality.To analyze the critical physical properties in LPBF, the following assumptions should be made to simplify the coupling physical phenomena.
1.The fluid in the molten pool is assumed to be Newtonian and laminar flow; 2. The mushy zone is assumed to be isotropic during the solid-liquid phase change process; 3. The distribution of powder particles is assumed to be fixed in space; 4. The effect of the surrounding argon's flow on the molten pool is neglected; 5.The influence of ambient gas on laser energy absorption is neglected.

Governing Equations
The mass conservation equation in LPBF is as follows [30]: where  ⃗ is the velocity of liquid metal; ρ is metal density.
For the non-uniform density distribution in the molten pool caused by thermal expansion, the buoyancy effects should be considered.Thus, the momentum conservation equation can be described as follows [31]: To analyze the critical physical properties in LPBF, the following assumptions should be made to simplify the coupling physical phenomena.

1.
The fluid in the molten pool is assumed to be Newtonian and laminar flow; 2.
The mushy zone is assumed to be isotropic during the solid-liquid phase change process; 3.
The distribution of powder particles is assumed to be fixed in space; 4.
The effect of the surrounding argon's flow on the molten pool is neglected; 5.
The influence of ambient gas on laser energy absorption is neglected.

Governing Equations
The mass conservation equation in LPBF is as follows [30]: where → u is the velocity of liquid metal; ρ is metal density.For the non-uniform density distribution in the molten pool caused by thermal expansion, the buoyancy effects should be considered.Thus, the momentum conservation equation can be described as follows [31]: where → F B is buoyancy described using the Boussinesq approximation → F m is used to characterize the difference of fluidity caused by the liquid-solid transition, which is calculated as force term caused by the solidification and can be described as D = CF 2 s / (1 − F s ) 3 + B ; F s is the local solid fraction in the computational domain; C is a constant that depends on the mushy zone microstructure; B is a small constant that is used to avoid D being infinity when Fs = 1.
The energy conservation equation is given as follows: where k is the thermal conductivity; q is the heat source; E is the specific internal energy which can be described as: where c is the specific heat; L is the latent heat of melting.
To track the free surface of molten pool, the volume of fluid (VOF) method is used, which can be described as: where γ is the volume fraction of liquid metal material in the computation element.

Initial and Boundary Conditions
In LPBF, the thermodynamic behavior at the boundary includes a laser heat source, heat dissipation by evaporation, heat convection, and radiation with the surrounding environment.So the thermal boundary condition and initial condition can be defined as follows [32]: T(x, y, z, 0) = T 0 (7) where T 0 is the ambient temperature; n is the vector along the normal direction of the top surface of the powder bed; q(x, y, z, t) represents the laser energy distribution; ε r is the surface radiation coefficient; σ s is the Stefan-Boltzmann constant (5.67 × 10 −8 W/(m 2 •K 4 )); h c is the natural convection heat transfer coefficient; q ev is the heat taken away as the molten pool evaporation.
Based on the characteristics of laser emission, the laser heat source is assumed to be Gaussian distribution [33,34].The heat flux of the laser decays exponentially from the center outward, and it is given as: where q is the heat flux; r is the radial distance from the laser center; θ is the laser radius where q = q 0 e −2 ; q 0 is the laser beam intensity at r = 0. Thus, q 0 can be written as: where ξ is the absorptivity of the powder bed; P is the laser power.Thus, the heat flux can be rewritten as: For high laser energy density, the molten pool temperature will exceed the evaporation temperature.The evaporative heat loss can be written as [35]: where R is the gas constant; P 0 is atmospheric pressure; ∆H* is the effective evaporation enthalpy; M is the molar mass; T v is the evaporation temperature.
As the evaporation of the molten pool under high laser energy density, the recoil pressure is considered as follows [31]: where L e is the latent heat of evaporation; T v is the evaporation temperature; P 0 is the atmospheric pressure.At the molten pool surface, the pressure balance can be described as [36]: where n is the normal direction of the molten pool surface; P is the pressure at the molten pool surface; µ is the dynamic viscosity; → v n is the normal velocity; R 1 and R 2 are the principal radii of molten pool surface curvature.The temperature-dependent surface tension coefficient can be described as σ = σ 0 − ∂σ/∂T T − T re f .
For the temperature dependence of the surface tension coefficient, the Marangoni effect appears.The corresponding fluid velocity boundary condition at the molten pool surface can be described as: where r is the tangential direction of the molten pool surface; → v t is the tangential velocity.

Simulation Cases and Computation
The simulation was carried out using the Flow3D V11.1 commercial software to simulate the melting-solidification process and thermal behavior.The simulation domain size is 1000 µm × 400 µm × 250 µm, and the powder layer thickness is 45 µm.The hexahedral mesh elements were used.The mesh size was set to be 5 µm, which is sufficient to describe the non-uniform powder distribution.The corresponding material properties of SS316L are shown in Table 1, and the temperature-dependent parameters are shown in Figure 2.

Molten Pool Morphology Analysis
When the molten pool temperature exceeded the evaporation temperature, the recoil pressure induced by the escaping vapor of the liquid metal was larger than its surface tension [41].Thus, the keyhole formed as shown in Figure 3.At the center of molten pool, the surface temperature was higher than that at the bottom.Thus, the buoyancy appeared due to thermal expansion, resulting in natural convection flowing from the bottom up.The temperature of the molten pool center was also higher than that at the edge.Thus, the surface tension at the molten pool center was smaller than that at the edge for the negative change rate of surface tension.Under the effects of the surface tension, the liquid metal at the center was drawn to the edge [42].The corresponding flow pattern is shown in Figure 3.This outward flow can promote the bonding of the adjacent scan track after solidification.The asymmetrical circulation in the molten pool was caused by the inhomogeneity of powder distribution, which also created the rough contour line of the scan track after the molten pool solidification.

Temperature Distribution Analysis
Linear energy density (LED) is commonly used to analyze the temperature distribution in LPBF [43]. Figure 4 shows the temperature distribution on the substrate along the scanning direction under the LED of 240 J/m, 400 J/m, 600 J/m, and 800 J/m respectively, which was obtained at the same laser power of 120 W. It can be found that when the molten pool temperature exceeded the evaporation temperature, the variation of LED had little effect on the temperature near the laser focus.This is because the input energy was taken away by evaporation [41], but the keyhole width increased significantly with the

Molten Pool Morphology Analysis
When the molten pool temperature exceeded the evaporation temperature, the recoil pressure induced by the escaping vapor of the liquid metal was larger than its surface tension [41].Thus, the keyhole formed as shown in Figure 3.

Molten Pool Morphology Analysis
When the molten pool temperature exceeded the evaporation temperature, the recoil pressure induced by the escaping vapor of the liquid metal was larger than its surface tension [41].Thus, the keyhole formed as shown in Figure 3.At the center of molten pool, the surface temperature was higher than that at the bottom.Thus, the buoyancy appeared due to thermal expansion, resulting in natural convection flowing from the bottom up.The temperature of the molten pool center was also higher than that at the edge.Thus, the surface tension at the molten pool center was smaller than that at the edge for the negative change rate of surface tension.Under the effects of the surface tension, the liquid metal at the center was drawn to the edge [42].The corresponding flow pattern is shown in Figure 3.This outward flow can promote the bonding of the adjacent scan track after solidification.The asymmetrical circulation in the molten pool was caused by the inhomogeneity of powder distribution, which also created the rough contour line of the scan track after the molten pool solidification.

Temperature Distribution Analysis
Linear energy density (LED) is commonly used to analyze the temperature distribution in LPBF [43]. Figure 4 shows the temperature distribution on the substrate along the scanning direction under the LED of 240 J/m, 400 J/m, 600 J/m, and 800 J/m respectively, which was obtained at the same laser power of 120 W. It can be found that when the molten pool temperature exceeded the evaporation temperature, the variation of LED had little effect on the temperature near the laser focus.This is because the input energy was taken away by evaporation [41], but the keyhole width increased significantly with the At the center of molten pool, the surface temperature was higher than that at the bottom.Thus, the buoyancy appeared due to thermal expansion, resulting in natural convection flowing from the bottom up.The temperature of the molten pool center was also higher than that at the edge.Thus, the surface tension at the molten pool center was smaller than that at the edge for the negative change rate of surface tension.Under the effects of the surface tension, the liquid metal at the center was drawn to the edge [42].The corresponding flow pattern is shown in Figure 3.This outward flow can promote the bonding of the adjacent scan track after solidification.The asymmetrical circulation in the molten pool was caused by the inhomogeneity of powder distribution, which also created the rough contour line of the scan track after the molten pool solidification.

Temperature Distribution Analysis
Linear energy density (LED) is commonly used to analyze the temperature distribution in LPBF [43]. Figure 4 shows the temperature distribution on the substrate along the scanning direction under the LED of 240 J/m, 400 J/m, 600 J/m, and 800 J/m respectively, which was obtained at the same laser power of 120 W. It can be found that when the molten pool temperature exceeded the evaporation temperature, the variation of LED had little effect on the temperature near the laser focus.This is because the input energy was taken away by evaporation [41], but the keyhole width increased significantly with the increase of LED for the evaporation enhancement.When the LED decreased to 240 J/m, the recoil pressure on the molten pool surface (Knudsen layer) was wakened, and the keyhole depth decreased significantly, which made it so the laser energy could not be sufficiently transferred to the substrate.Thus, the substrate temperature was lower than the liquidus temperature, which means the previously processed layer and the currently processed one could not be effectively fused in the product manufacturing process, resulting in poor densification [8].If the LED was further decreased, the densification would drop sharply due to the incomplete melting of the powders [44].
increase of LED for the evaporation enhancement.When the LED decreased to 240 J/ the recoil pressure on the molten pool surface (Knudsen layer) was wakened, and keyhole depth decreased significantly, which made it so the laser energy could not sufficiently transferred to the substrate.Thus, the substrate temperature was lower th the liquidus temperature, which means the previously processed layer and the curren processed one could not be effectively fused in the product manufacturing process, resu ing in poor densification [8].If the LED was further decreased, the densification wou drop sharply due to the incomplete melting of the powders [44].To analyze the substrate temperature at laser focus, a temperature sampling pro was set on the substrate, which moved with the laser focus synchronously.To well obser the temperature on the substrate, the keyhole formation on the substrate should avoided.Thus, the laser power = 120 W and scanning speed = 0.45 m/s were used in t simulation.The temperature at the probe is shown in Figure 5.It could be found that wh the laser turned on, the temperature rose quickly.At t = 0.15 ms, the probe temperatu reached the liquidus temperature, and the molten pool formed.For the energy absorpti of powder melting, the probe temperature of the molten pool rose slowly until the mol pool formed sufficiently.Then, the probe temperature fluctuated around 2635 K (dash line) due to the non-uniform powder distribution.Moreover, a fixed temperature sampling probe was set up on the substrate to analy the multi-track scanning process.In this simulation, the laser power = 120 scanni speed = 0.45 m/s, and hatch space = 0.12 mm.The corresponding scanning strategy shown in Figure 6, and the evolution of temperature is shown in Figure 7.It can be se that the sampling temperature decreased periodically.In each period, the temperatu rose rapidly when the laser focus was close to the sampling probe, and the temperatu decreased relatively slowly when the laser focus was far away from the sampling pro To analyze the substrate temperature at laser focus, a temperature sampling probe was set on the substrate, which moved with the laser focus synchronously.To well observe the temperature on the substrate, the keyhole formation on the substrate should be avoided.Thus, the laser power = 120 W and scanning speed = 0.45 m/s were used in this simulation.The temperature at the probe is shown in Figure 5.It could be found that when the laser turned on, the temperature rose quickly.At t = 0.15 ms, the probe temperature reached the liquidus temperature, and the molten pool formed.For the energy absorption of powder melting, the probe temperature of the molten pool rose slowly until the molten pool formed sufficiently.Then, the probe temperature fluctuated around 2635 K (dashed line) due to the non-uniform powder distribution.
the liquidus temperature, which means the previously processed layer and the cu processed one could not be effectively fused in the product manufacturing process ing in poor densification [8].If the LED was further decreased, the densification drop sharply due to the incomplete melting of the powders [44].To analyze the substrate temperature at laser focus, a temperature samplin was set on the substrate, which moved with the laser focus synchronously.To well the temperature on the substrate, the keyhole formation on the substrate sho avoided.Thus, the laser power = 120 W and scanning speed = 0.45 m/s were used simulation.The temperature at the probe is shown in Figure 5.It could be found th the laser turned on, the temperature rose quickly.At t = 0.15 ms, the probe temp reached the liquidus temperature, and the molten pool formed.For the energy abs of powder melting, the probe temperature of the molten pool rose slowly until the pool formed sufficiently.Then, the probe temperature fluctuated around 2635 K line) due to the non-uniform powder distribution.Moreover, a fixed temperature sampling probe was set up on the substrate to the multi-track scanning process.In this simulation, the laser power = 120 W, sc speed = 0.45 m/s, and hatch space = 0.12 mm.The corresponding scanning stra shown in Figure 6, and the evolution of temperature is shown in Figure 7.It can that the sampling temperature decreased periodically.In each period, the temp rose rapidly when the laser focus was close to the sampling probe, and the temp decreased relatively slowly when the laser focus was far away from the sampling Moreover, a fixed temperature sampling probe was set up on the substrate to analyze the multi-track scanning process.In this simulation, the laser power = 120 W, scanning speed = 0.45 m/s, and hatch space = 0.12 mm.The corresponding scanning strategy is shown in Figure 6, and the evolution of temperature is shown in Figure 7.It can be seen that the sampling temperature decreased periodically.In each period, the temperature rose rapidly when the laser focus was close to the sampling probe, and the temperature decreased relatively slowly when the laser focus was far away from the sampling probe.The peak temperature in each period was extracted and fitted as the dashed line in Figure 7, which can be described using the following formula: where t(ms) is laser scanning time.Based on this equation, it can be predicted that the temperature can be reduced to 320 K when the scanning time is 556 ms.
The peak temperature in each period was extracted and fitted as the dashed line in 7, which can be described using the following formula: where t(ms) is laser scanning time.Based on this equation, it can be predicted temperature can be reduced to 320 K when the scanning time is 556 ms.

Keyhole Depth Analysis
The keyhole has great effects on the pore formation, which may degrade the mance of the final parts [45].To simulate the keyhole formation process, differe scanning speeds and scanning powers were analyzed.Figure 8 shows the molt surface distribution for four different scanning speeds (0.15, 0.2, 0.3, and 0.4 m/s same laser power of 120 W, respectively.Figure 9 shows the molten pool surface d tion for five different laser powers (90, 105, 120, 135, and 150 W) for the same s speed of 0.2 m/s, respectively.Then, the keyhole depth was obtained, as shown in 10 and Figure 11, respectively.It could be found that with the decrease of scannin and the increase of laser power, the gradient of the keyhole depth also increas indicates that the keyhole depth can be increased more effectively by increasing ergy input after the keyhole has formed.Moreover, the remelting of the previo can increase the densification and wetting properties [44].The peak temperature in each period was extracted and fitted as the dashed lin 7, which can be described using the following formula: where t(ms) is laser scanning time.Based on this equation, it can be predicte temperature can be reduced to 320 K when the scanning time is 556 ms.

Keyhole Depth Analysis
The keyhole has great effects on the pore formation, which may degrade t mance of the final parts [45].To simulate the keyhole formation process, diff scanning speeds and scanning powers were analyzed.Figure 8 shows the m surface distribution for four different scanning speeds (0.15, 0.2, 0.3, and 0.4 m same laser power of 120 W, respectively.Figure 9 shows the molten pool surfac tion for five different laser powers (90, 105, 120, 135, and 150 W) for the same speed of 0.2 m/s, respectively.Then, the keyhole depth was obtained, as shown 10 and Figure 11, respectively.It could be found that with the decrease of scann and the increase of laser power, the gradient of the keyhole depth also incre indicates that the keyhole depth can be increased more effectively by increasi ergy input after the keyhole has formed.Moreover, the remelting of the prev can increase the densification and wetting properties [44].

Keyhole Depth Analysis
The keyhole has great effects on the pore formation, which may degrade the performance of the final parts [45].To simulate the keyhole formation process, different laser scanning speeds and scanning powers were analyzed.Figure 8 shows the molten pool surface distribution for four different scanning speeds (0.15, 0.2, 0.3, and 0.4 m/s) for the same laser power of 120 W, respectively.Figure 9 shows the molten pool surface distribution for five different laser powers (90, 105, 120, 135, and 150 W) for the same scanning speed of 0.2 m/s, respectively.Then, the keyhole depth was obtained, as shown in Figures 10 and 11, respectively.It could be found that with the decrease of scanning speed and the increase of laser power, the gradient of the keyhole depth also increased.This indicates that the keyhole depth can be increased more effectively by increasing the energy input after the keyhole has formed.Moreover, the remelting of the previous track can increase the densification and wetting properties [44].
The peak temperature in each period was extracted and fitted as the dashed line in Figure 7, which can be described using the following formula: where t(ms) is laser scanning time.Based on this equation, it can be predicted that the temperature can be reduced to 320 K when the scanning time is 556 ms.

Keyhole Depth Analysis
The keyhole has great effects on the pore formation, which may degrade the performance of the final parts [45].To simulate the keyhole formation process, different laser scanning speeds and scanning powers were analyzed.Figure 8 shows the molten pool surface distribution for four different scanning speeds (0.15, 0.2, 0.3, and 0.4 m/s) for the same laser power of 120 W, respectively.Figure 9 shows the molten pool surface distribution for five different laser powers (90, 105, 120, 135, and 150 W) for the same scanning speed of 0.2 m/s, respectively.Then, the keyhole depth was obtained, as shown in Figure 10 and Figure 11, respectively.It could be found that with the decrease of scanning speed and the increase of laser power, the gradient of the keyhole depth also increased.This indicates that the keyhole depth can be increased more effectively by increasing the energy input after the keyhole has formed.Moreover, the remelting of the previous track can increase the densification and wetting properties [44].To analyze the multiple-track interaction, the multiple-track manufacturing process was simulated with laser power = 120 W, scanning speed = 0.2 m/s, and hatch space = 0.12 mm.The molten pool surface distribution of the second and third tracks is shown in Fig ure 12.It shows that the two molten pools had the same morphology, which means tha the residual heat from the previous track had little effect on the keyhole depth of the cur rent one.To analyze the multiple-track interaction, the multiple-track manufacturing process was simulated with laser power = 120 W, scanning speed = 0.2 m/s, and hatch space = 0.12 mm.The molten pool surface distribution of the second and third tracks is shown in Figure 12.It shows that the two molten pools had the same morphology, which means that the residual heat from the previous track had little effect on the keyhole depth of the current one.To analyze the multiple-track interaction, the multiple-track manufacturing process was simulated with laser power = 120 W, scanning speed = 0.2 m/s, and hatch space = 0.12 mm.The molten pool surface distribution of the second and third tracks is shown in Figure 12.It shows that the two molten pools had the same morphology, which means that the residual heat from the previous track had little effect on the keyhole depth of the current one.To analyze the multiple-track interaction, the multiple-track manufacturing process was simulated with laser power = 120 W, scanning speed = 0.2 m/s, and hatch space = 0.12 mm.The molten pool surface distribution of the second and third tracks is shown in Figure 12.It shows that the two molten pools had the same morphology, which means that the residual heat from the previous track had little effect on the keyhole depth of the current one.To analyze the multiple-track interaction, the multiple-track manufacturing process was simulated with laser power = 120 W, scanning speed = 0.2 m/s, and hatch space = 0.12 mm.The molten pool surface distribution of the second and third tracks is shown in Figure 12.It shows that the two molten pools had the same morphology, which means that the residual heat from the previous track had little effect on the keyhole depth of the current one.

Conclusions
In this study, the multi-physics simulation model was developed to analyze the melting-solidification process of SS316L in the laser powder bed fusion.Using this model, the temperature distribution and molten pool morphology were analyzed by changing the process parameters.Moreover, the effects of powder distribution on workpiece quality were analyzed.Some valuable conclusions were drawn, as follows: 1.
The convection flow in the molten pool can effectively widen the molten pool width and promote strong bonding between adjacent scan tracks.Therefore, the hatch space can be enlarged by increasing the laser power or decreasing the scanning speed to enhance the convection flow behavior.

2.
When the LED decreases to 240 J/m, the keyhole depth becomes too small to fuse the previously processed layer with the currently processed one, potentially leading to degradation in part densification.To ensure better densification of the final parts, it is suggested that the LED be set to over 400 J/m when the layer thickness is 45 µm.

3.
The keyhole depth can be enlarged more effectively by further increasing the energy input after the keyhole is formed, such as by increasing the laser power or decreasing the scanning speed.

4.
If the hatch space is too large or the powder bed is sparsely distributed, internal void defects may form, significantly affecting workpiece quality.To prevent these defects, it is suggested that the hatch space be narrower than the single-track width.
Furthermore, the comprehensive liquid-gas phase change model was not considered for the obstacle of integration with the current simulation model.This may result in the deviation of the temperature field between our numerical simulation and experiments when reaching evaporation temperature.In the follow-up research, a sophisticated simulation model including the liquid-gas phase change should be developed to comprehensively analyze the dynamics of the temperature field after the appearance of evaporation.

Figure 1 .
Figure 1.Schematic diagram of the LPBF process.The gray represents the unprocessed steel; the blue represents the molten pool; and the baby blue represents the molten pool after solidification.

Figure 1 .
Figure 1.Schematic diagram of the LPBF process.The gray represents the unprocessed steel; the blue represents the molten pool; and the baby blue represents the molten pool after solidification.

Figure 3 .
Figure 3. Flow behavior in the molten pool with different cross sections (laser power = 120 W, scanning speed = 0.2 m/s).The black arrows show the fluid flow; the arrow length shows the velocity magnitude; and the color scale shows the temperature field in Kelvin.

Figure 3 .
Figure 3. Flow behavior in the molten pool with different cross sections (laser power = 120 W, scanning speed = 0.2 m/s).The black arrows show the fluid flow; the arrow length shows the velocity magnitude; and the color scale shows the temperature field in Kelvin.

Figure 3 .
Figure 3. Flow behavior in the molten pool with different cross sections (laser power = 120 W, scanning speed = 0.2 m/s).The black arrows show the fluid flow; the arrow length shows the velocity magnitude; and the color scale shows the temperature field in Kelvin.

Figure 7 .
Figure 7. Evolution of substrate temperature at the fixed sampling probe.

Figure 8 .
Figure 8. Molten pool surface distribution under different scanning speeds.

Figure 7 .
Figure 7. Evolution of substrate temperature at the fixed sampling probe.

Figure 8 .
Figure 8. Molten pool surface distribution under different scanning speeds.

Figure 7 .
Figure 7. Evolution of substrate temperature at the fixed sampling probe.

Figure 7 .
Figure 7. Evolution of substrate temperature at the fixed sampling probe.

Figure 8 .
Figure 8. Molten pool surface distribution under different scanning speeds.

Figure 12 .
Figure 12.Molten pool surface distribution during multiple-track fabrication.

Figure 12 .
Figure 12.Molten pool surface distribution during multiple-track fabrication.

Figure 12 .
Figure 12.Molten pool surface distribution during multiple-track fabrication.

Figure 12 .
Figure 12.Molten pool surface distribution during multiple-track fabrication.

Figure 12 .
Figure 12.Molten pool surface distribution during multiple-track fabrication.

Table 1 .
Thermal and physical properties of SS316L.